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Abstract 

In the very early stages of LHC running, uncertainties in detector performance will lead to large 
ambiguities in jet, electron and photon energy measurements, along with inferred missing transverse 
energy EJp^^^. However, muon detection should be quite straightforward, with the added benefit 
that muons can be reliably detected down to transverse energies of order 5 GeV. Supersymmetry 
discovery through multimuon channels has been extensively explored in the literature, but always 
relying on hard EJp"^^^ cuts. Here, we quantify signal and background rates for same-sign (SS) 
dimuon and multimuon production at the LHC without any EJp^^'^ cuts. The LHC, operating at 
^/s = 10 TeV, should be able to discover a signal over expected background consistent with gluino 
pair production for mg ~ 450 (550) GeV in the SS dimuon plus > 4 jets state with just 0.1 (0.2) 
fb^^ of integrated luminosity. 

PACS numbers: ll.lO.Kk, 14.80.-j, 04.50.+h 
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With the recent circulation of proton beams around the entire CERN Large Hadron 
CoUider (LHC) ring, the era of LHC physics has begun. Meaningful data is now expected 
starting in fall 2009, when LHC will likely start up with pp collisions at ^/s ~ 10 TeV. 

In the LHC ramping up process, it will be essential to observe many familiar Standard 
Model (SM) processes- multi-jet production as predicted by QCD, W and Z production 
as predicted by the electroweak theory, ti production, vector boson pair production- all at 
their expected rates, and with distributions and mass peaks at ,evio„.y.ea..edvaluesQ. 
Conventional wisdom holds that once confidence in the Atlas and CMS detectors has been 
established, then the search for physics beyond the Standard Model will begin. In this letter, 
we explore the possibility of searching for new physics in parallel with the calibration phase. 
We will show that even with relatively poor knowledge of the detector, new physics searches 
may still be possible, at least in the case of weak scale supersymmetry. 

Weak scale supersymmetry- wherein each particle state of the Standard Model has a TeV- 
scale superpartner differing by 1/2/?, units of spin- is perhaps the most motivated new physics 
theory [2]. Theories with supersymmetry (SUSY) broken at the weak scale actually enjoy 
indirect experimental support in that the measured values of the three SM gauge couplings 
at energy scale Q ~ Mz, when extrapolated to very high energies under renormalization 
group (RG) evolution, meet at a point as predicted by grand unified theories (GUTs) under 
Minimal Supersymmetric Standard Model (MSSM) evolution (while they miss badly under 
SM evolution) [31] . SUSY theories also predict a SM-like Higgs boson h with mass below 
~ 135 GeV- a scenario which is consistent with global analyses of precision electroweak 
measurements jj] . 

While the idea of supersymmetry is theoretically very appealing, the mechanism behind 
SUSY breaking is a complete mystery. One very elegant SUSY breaking mechanism occurs 
in local SUSY- or supergravity (SUGRA)- theories. It is possible to embed the SM into 
a supergravity theory, and then set up a hidden sector which serves as an arena for SUSY 
breaking via the super-Higgs mechanism. The SUSY breaking is communicated from the 
hidden sector to the visible sector via Planck-scale suppressed operators, and a judicious 
choice of parameters leads to weak scale soft SUSY breaking (SSB) parameters, exactly 
as needed by gauge coupling evolution, and which serve to stabilize the weak scale-GUT 
scale energy hierarchy without too much fine-tuning. The simplest such model, the minimal 
supergravity or mSUGRA model, thus posits a common (universal) mass mo for all SSB 
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scalar masses, a common SSB gaugino mass mi/2, and common trilinear SSB terms A^. 
(Here, the gaugino is the spin-| superpartner of the gauge bosons.) Motivated by gauge 
couphng unification, these common masses are assumed vahd at the GUT scale Mqut — 
2 X 10^^ GeV. All weak-scale Lagrangian parameters can be calculated in terms of this 
parameter set using the power of the RG equations. Thus, all physical superpartner masses 
and mixings may be calculated in terms of the parameter set 

mo, mi/2, v4o, tan (3 , sign{^) , (1) 

wherein tan (3 is the ratio of the two Higgs field vacuum expectation values (vevs) needed 
for electroweak symmetry breaking, and is a quadratic superpotential term. While many 
other well-motivated SUSY models exist, the mSUGRA model has emerged as a sort of 
paradigm choice for exploring basic SUSY phenomena expected at collider experiments. 

The strongly interacting sparticles- the gluinos g and squarks q- often end up with the 
largest of all the sparticle masses due to the infiuence of the strong interactions on their RG 
mass evolution. Sparticles such as charginos, neutralinos and sleptons are frequently much 
lighter. The strongly interacting g and q- produced through QCD interactions- usually have 
the largest production cross sections. Once the g and gs are produced, they decay through a 
cascade of possibly several stages until the state with the lightest SUSY particle- or LSP- is 
reachedj^. The LSP in mSUGRA usually turns out to be the lightest neutralino, Zi, which 
if i?-parity is conserved, is absolutely stable and serves as a good candidate for cold dark 
matter (CDM) in the universelsl. 

The classic signature for g and q production at hadron colliders consists of events con- 
taining jets plus large missing transverse energy E^'^'^^, wherein the E^^^^ arises due to the 
Zis completely escaping the detector, much as neutrinos do. This signature channel should 
serve sparticle-hunters well once the detectors are fully calibrated so that SM backgrounds 
for jets+i?^**'^ events are well-understood [3]. Experience with similar jets+E'^*'''^ searches at 
the Fermilab Tevatron suggest that it may well take some time to fully understand detector 
performance, so that i?™'^* can be reliably measured. For this reason, several of us recently 
proposed that early searches for SUSY matter at the LHC may be better served by looking 
for events containing multiple (2,3,4,...), high transverse momentum {pt > 20 GeV) isolated 
leptons (es and/or /is) along with jets, instead of i?™'^*+jets events^]. Requiring high lepton 
multiplicity rejects SM background at a large rate, while maintaining much of the expected 
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signal, since isolated leptons are expected to be produced frequently in the sparticle cascade 
decays [s]. 

Since publication of Ref. 81], it has been pointed out that reliable electron identification 
may also be a major issue during the early phase of LHC running. If so, this could jeopardize 
the results of Ref. Q], which summed over both muons and electrons in order to establish 
the multi-lepton signal and background rates. In addition, the SM background calculation 
of Ref. 8| included only 2^2 processes that were pre-programmed into Isajet. However, 
various SM 2 — n BG processes potentially may be larger than the lowest order processes 



considered in Ref. 

In this letter, we show i). that it is sufficient to focus only on isolated multimuon plus 
jets events during the earliest SUSY searches at LHC The lack of electron channels can 
be partially compensated for by the lower pt values which are allowed for isolated muon 
searches. Secondly, ii). we evaluate a variety of additional 2 — > n background processes 
beyond those presented in Ref. [81], thus putting our results on a more firm foundation. 
Thirdly, we re-evaluate all signal and background channels for the anticipated start-up energy 
of \/s = 10 TeV, instead of design energy ^/s = 14 TeV. Finally, in), we scan over a wide 
swath of mSUGRA model parameter space, and present the LHC reach plot on the m(squark) 
vs. m(gluino) plane for various low levels of integrated luminosity. In the same-sign dimuon 
plus jets channel, some reach is possible even for integrated luminosities as low as 0.1 fb~^, 
where squark and gluino masses up to ~ 450 GeV may be probed. 

There are several advantages to a SUSY search via multimuon plus jets events. 

• Reliable electron identification may be difficult in the early stages of LHC running. 
Electrons will need to be readily distinguishable from QCD jets and also from high 
photon production. As an example, a jet with a single soft charged pion plus several 
7r°s can give a track pointing to a mainly electromagnetic calorimeter deposition, which 
may well fake an electron signal. 

• Muon identification should be straightforward even in the very early stages of LHC 
running 10(]. In fact, cosmic ray muons have already been seen at both Atlas and 



CMS. Muons with ~ 5 GeV should readily penetrate the electromagnetic calorime- 
ter (ECAL) and hadronic calorimeter (HCAL), yielding easily-seen tracks in the muon 
chambers. Since muons are so heavy, they produce minimal bremsstrahlung and show- 
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ering in the ECAL or HCAL. 

• Muons can be readily identified at pt values much lower than electrons. Reliable e 
tagging typically needs pri^) ~ 20 GeV, while PtI/u) ~ 5 GeV is sufficient for muon 
identification. Thus, the lower pt muons emerging from cascade decays will be easily 
detected, while this is not so for electrons. 

• Superparticle cascade decays tend to be rich in 5s and rs. While b ciiv^ decay yields 
mainly non-isolated muons, r — > v^-^v^ decay leads to rather soft, but isolated, muon 
production. The rather low Prif^) requirements allows one to detect muons from r 
decay, while es from tau decay are often too soft to reliably identify. 

The large rate for b and r production in sparticle cascade decay events has three 



sources llj: 1. the large b and r Yukawa couplings, especially at large tan/3 values, en- 
hance chargino and neutralino branching fractions into b and r states, 2. third generation 
sparticle masses are often much lighter than their ffist / second generation counterparts due 
to Yukawa coupling effects pushing the third generation SSB masses to low values, and 
also due to large mixing effects, which are proportional to the corresponding fermion mass. 
This latter effect also enhances sparticle decay rates into third generation fermions. 3. Higgs 
bosons, especially h, can be produced at large rates in sparticle cascade decays. For instance, 
if the decay Z2 Zih is kinematically allowed, this usually dominates the Z2 branching 
fraction. Since h and the other Higgs subsequently decay dominantly into third generation 
fermions, one gets enhanced b and r production from cascade decays of sparticles into Higgs 
bosons. 

The search for multi-muon events has been proposed much earlier with regards to the 



search for fourth generation quarks[l2j, which also decay via a cascade to the lightest flavor 
states. Multimuon detection has been proposed in the old idea of an "iron ball detector", 
wherein the interaction region^ completely surrounded by iron absorber, and one only 
detects the penetrating muons[13]. LHC detectors are vastly more complex than the iron 
ball detector. But in the very early stages of running, wherein calorimeter and other detector 
response is not well understood, their initial performance may approximate the iron-ball idea. 

We adopt the Isajet 7.78 program for sparticle mass calculations and simulation of signal 
events at the LHC [14]. A toy detector simulation is employed with calorimeter cell size 
Ar] X A0 = 0.05 x 0.05 and —5<r]<5. (here, 77 = —log tan | is pseudorapidity and 



is angle transverse to beamline). The HCAL (hadronic calorimetry) energy resolution is 
taken to be 80%//E + 3% for \r]\ < 2.6 and FCAL (forward calorimetry) is 100%//E + 5% 
for \ri\ > 2.6. The ECAL (electromagnetic calorimetry) energy resolution is assumed to 
be 'i%j\fE + 0.5%. We use the Isajet jl^ jet finding algorithm with jet cone size R = 
+ A02 = 0.4 and require that ET{jet) > 50 GeV and \i]{jet)\ < 3.0. Muons are 
considered isolated if they have prif^) > 5 GeV and |?7^| < 2 with visible activity within 
a cone of AR < 0.2 of TiE^^^'^ < 5 GeV. The isolation criterion helps reduce multi-lepton 
backgrounds from heavy quark (cc and hb) production. 

For our initial analysis, we adopt the well-studied SPSla' benchmark point [l^, which 
occurs in the minimal supergravity (mSUGRA) model with parameters mo = 70 GeV, 
mi/2 = 250 GeV, Aq = -300 GeV, tan/5 = 10, /i > and mt = 172.6 GeV. The SPSla' 
point leads to a spectrum with rrig = 608 GeV, while squark masses tend to be in the 550 
GeV range. The gluinos and squarks then cascade decay via a multitude of modes leading 
to events with high jet, 6-jet, isolated lepton and tau lepton multiplicity. 

Since the gluino and squark cascade decay events will be rich in jet activity, we first 
require events with > 4 jets, with Exijl, j2, jS, ji) > 100, 50,50,50 GeV. We also require 
sphericity (restricted to the transverse plane) St > 0.2 to reject QCD-like events at little 
cost to signal. We do not apply the traditional cut on missing transverse energy, since at 
this stage we are working towards early SUSY discovery, when EJp^'^'^ may not yet be well- 
established. In Fig. [T] we show the muon px distribution from point SPSla' along with 
dominant SM BGs (see discussion below) in same-sign (SS) dimuon plus > 4 jet events 
(each muon in a SS-dimuon event will have an entry in the plot). The signal muons tend 
to populate the 5 — 40 GeV regime, and so should be easily measured by the bend of their 
tracks in the detector magnetic field. The BG muons come dominantly from tt production, 
and have a hard component (from W decay) and a soft component (from rare h and c decays 
to isolated muons). The soft component exceeds signal in the 5-10 GeV range. Hence, we 
require prif^) > 10 GeV in our multi-muon signal events, which eliminates much of the BG 
from ti production. 

Next, we plot the multiplicity of muons in the SUSY cascade decay events. The results 
are shown for point SPSla' in Fig. [2]for pp collisions at ^/s = 10 TeV. We also plot a variety 
of SM backgrounds. The dominant backgrounds for the dimuon signal were calculated using 
AlpGen/Pythia and AlpGen's matching algorithm (MLM scheme 16|]) to include multiple jet 
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FIG. 1: pr(A*) distribution from the SPSla' mSUGRA study, along with SM BG, for SS dimuon 
plus > 4 jet events at LHC with = 10 TeV. 

emission. In particular, the ti channel includes tt + 0, 1, 2, 3 jets, Z+jet includes Z + 0, 1, 2, 3 
jets, the ti + Z channel includes ti + Z + 0,1,2 jets and bb + Z includes bb + Z + 0,1,2 
jets (in all cases the full matrix element 'j*,Z* —>■ l^l^ was used). The presence of hard 
additional jets increases the BG quite a bit from our earlier estimates using just the Isajet 
parton shower.^ 

In addition, we have calculated using MadGraph/Pythia jl?!. fisl] a variety of exact 2 ^ n 
processes: titi, tibb, bbbb, tiV, tiVV, VV, W + jet, bb, QCD dijets, VVV and WW 
production, where V = or The summed BG histogram along with component 
contributions are also shown in Fig. [2l We adopt a renormalization/factorization scale 
choice Q = {s is the parton-parton CM frame squared energy) which brings our 

background (BG) cross sections into close accord with NLO QCD results.^ 



^ Our BG from bb production comes from Alpgen/Pythia, but with just LO bb production along with jets 
from the parton shower. We expect the BG from 66+jets production to be sub-dominant because we 
would need to obtain one isolated lepton from a b decay, and another from a c decay, while producing 
four hard jets at the same time, and that sort of event is extremely rare. This reaction, with exact 4-jet 
emission matrix elements, is extremely hard to generate with reliable statistics. 



^ The dominant BG for SS dimuon plus > 4-jet events comes from tt production. Using the MCFM code 19|, 
we find a^'-'{pp ^ tt){Q — nitop) — 255 pb, while a^^'~'{pp ^ tt){Q — rritop) — 347 pb. If instead we take 
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FIG. 2: Muon multiplicity cross sections expected from the SPSla' mSUGRA case study, along 
with SM background, for > 4 jet plus n-muon events (Pt(/^) ^ 10 GeV) at LHC with ^/s = 10 
TeV. 

At n(/i) = 0, signal is about three orders of magnitude below SM background. As we 
increase the isolated muon multiplicity, BG falls off faster than signal, and signal exceeds 
BG already at n{fj^) = 3, where signal is at the ~ 5 fb level. The dimuon signal can be 



2l|. 



broken up into opposite sign /i^/i" events (OS) and same sign events (SSljol. H. 

In order to suppress the large contribution from the Z peak to the OS events we required 
10 GeV< m{fi^fi^) < 75 GeV. The SS signal is due in part to the Majorana nature of the 
gluinos, in that a gluino is as likely to decay via g — > Wiud as via the charge conjugate mode. 

Thus, gg production is likely to lead to equal amounts of ++ and SS dileptons. Now, 

gq or qq production depends on the quark content of the colliding beams, and since LHC is 

a pp collider, we expect more ++ dileptons than dileptons from squark production. 

For OS dimuons and case study SPSla' , the OS signal is just slightly above OS BG, while 
SS signal well exceeds SS BG, and is at the ~ 10 fb level. As we move to higher and higher 
muon multiplicity, the signal rates diminish, although signal-to-background ratio steadily 
improves. For instance, at n{fi) = 3, signal is ~ 5 fb, while the summed SM background, 
arising mainly from ti and tiZ production, occurs at the ~ 0.1 fb level. Using these results. 



Q = \/j/6, then a^^{pp ^ ttX){Q = y/S/6) ~ 337 pb. 
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FIG. 3: Plot of a = pr(/i2)/"T'(/U^Ai^) for SS dimuon plus jets events at LHC with ^/s = 10 TeV, 
after cuts listed in text. 

we can now see that if case study SPSla' describes SUSY, then the first clear signal may 
emerge in the SS dimuon plus multi-jet channel, with corroborating signals in the OS and 
tri-muon channels. 

Another discriminating variable for SUSY events has been proposed by Randall and 



Tucker-Smith 2^], albeit applied to dijet events coming from squark pair production. They 
propose using a = pT{jet2) /m{jeti, jet2). Here we plot = pt(/^2)/'^(/^^/^^) in Fig. [3l 
and do find that the SUSY event shape is discriminated from the SM event shape. We do 
not apply an a{fi) cut at this time for very low luminosity studies, but merely note that this 
distribution will add additional confidence in any possible SS dimuon signal. 

As higher integrated luminosities are reached, trimuon and later four muon plus jet events 
should emerge at rates far above expected background. Of course, also as higher integrated 
luminosities are achieved, reliable electron ID should become available, and ultimately also 
reliable E^'^^^ measurements. Thus, the real utility of multi-muon plus jets events will be 
for a possible early discovery of SUSY, when muon ID is possible, but electron ID and EJp^^'^ 
resolution are still works in progress. 

In Fig. m we scan over ~ 200 choices of mo and mi/2 values for fixed Aq = 0, fi > 
and tan /5 = 45. We test to see if the SS dimuon plus jets signal is greater than a nominal 
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FIG. 4: Reach of the = 10 TeV LHC for mSUGRA models with Aq = {) and tan /3 = 45 via SS 
dimuon + > 4 jet events in the nig vs. rriq plane, for various integrated luminosity values. 



discovery threshold of 5a, and require at least five signal events as well, for various integrated 
luminosity choices: 0.1, 0.2 and 1 fb~^. We plot the results in the physical rrig vs. m^^ 
plane. The lower-right region gives a chargino mass less than 103.5 GeV, and so is already 
excluded by LEP2 new particle searches. The left side of the plot gives a fi slepton as the 
LSP, and is excluded by null searches for stable, charged relics from the Big Bang. For 
just 0.1 fb~^ of integrated luminosity, eleven points are accessible, with rrig ~ 480 GeV and 
~ 580 GeV. For 0.2 fb~^, rrig ~ 550 and ~ 700 GeV are being probed. The SS 
di-muon reach increases to nig ~ 650 GeV for 1 fb~^ of integrated luminosity. If we move to 
a lower tan (3 = 10 value, then the dimuon reach diminishes only slightly from that presented 
in the tan /5 = 45 case. 

Conclusions: In the early stages of LHC running, electron ID, ECAL and HCAL calibra- 
tion and EJp'^'^''^ resolution may all be works in progress. However, muon ID and momentum 
resolution, obtained from track bending in the magnetic field, should be quite reliable, and 
allow muon pj- measurement down to the ~ 5 — 10 GeV range. SS dimuon and, later, OS 
dimuon and > 3fi plus multi-jet signals, without any E^^^^ discrimination, should allow for 
good signal-to-background resolution for gluino masses up to about 550 GeV with just 0.2 
fb~^ of data. Thus, SS dimuon and multi-muon plus jets production offer excellent possi- 
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bilities for an early SUSY discovery at LHC, even if E^^^^ and electron ID are not initially 
well- understood. 
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